Mass distribution exponents for growing trees 



F. David 1 , P. Di Francesco 1 , E. Guitter 1 and T. Jonsson 2 

1 Service de Physique Theorique, CEA/DSM/SPhT 
Unite de recherche associee au CNRS 
CEA/Saclay 
91191 Gif sur Yvette Cedex, 
France 

2 Science Institute, 
University of Iceland, 
Dunhaga 3, 107 Reykjavik, 
Iceland 

f rancois . david@cea . f r 
philippe . di-francesco@cea.fr 
emmanuel.guitter@cea.fr 
thj onsOraunvis . hi . is 



Abstract 

We investigate the statistics of trees grown from some initial tree by attaching 
links to preexisting vertices, with attachment probabilities depending only on 
the valence of these vertices. We consider the asymptotic mass distribution 
that measures the repartition of the mass of large trees between their differ- 
ent subtrees. This distribution is shown to be a broad distribution and we 
derive explicit expressions for scaling exponents that characterize its behavior 
when one subtree is much smaller than the others. We show in particular the 
existence of various regimes with different values of these mass distribution 
exponents. Our results are corroborated by a number of exact solutions for 
particular solvable well as by numerical simulations. 



1. Introduction 

Random trees arise in many branches of science, ranging from the social sciences 
through biology, physics and computer science to pure mathematics. In physics, random 
trees often occur in the context of statistical mechanics or quantum field theory where 
the weight or probability of a tree is given by a local function, e.g. a Boltzmann factor 
with an energy which depends only on the valence of individual vertices but not on 
global features of the tree. A motivation for study of such trees comes for instance 
from results in the theory of random surfaces which behave like trees in some cases, see 
. Another motivation concerns the study of two-dimensional quantum gravity by use 
of so-called causal tessellations @, directly expressible in terms of random walks or of 
trees 0. 

Another large class of random trees arises in growth processes where a tree evolves in 
time by adding new vertices which attach themselves by a link to one of the vertices of a 
preexisting tree according to some stochastic rules, with local attachment probabilities. 
A possible realization is a tree-like molecule (branched polymer) growing in a solution 
with an abundance of monomers. In many important real world cases one can observe 
the structure of the growing tree but only guess the rules which govern the growth 
(social networks, the internet, citation networks etc.). For a review of this topic, see Pj. 

Most tree models in statistical mechanics with a local energy function fall into one 
universality class, called generic random trees, with susceptibility exponent 7 = 1/2, 
intrinsic Hausdorff dimension 2 and spectral dimension 4/3 ||. This universality class 
corresponds also to the mean field theory of branched polymers. If we look at a generic 
infinite tree, then with probability 1 there is a unique infinite non-backtracking path . 
A rooted generic tree can therefore be viewed as an infinite half-line with finite trees 
growing out of it and these finite outgrowths have a well-understood distribution. This 
means that if we sit at a vertex of a generic random tree, almost all the vertices of the 
tree are to be found in the direction of one of the links emanating from this vertex. 

There do not seem to be any local growth rules for a tree which produce generic 
trees. In this paper we study the structure of growing trees and focus on how they differ 
from generic trees. Taking binary planar trees with a root of valence one as an example 
we can ask what proportion of the vertices sits on the right hand side of the tree and 
what proportion on the left. For a generic tree almost all the vertices are on one side. We 
will see that with local growth rules, in the limit of infinite trees, we get a continuous 
distribution for the proportion of vertices sitting on each side. We will refer to this 
distribution as the mass distribution. This quantity is one of the simplest geometrical 
characterizations of the geometry of the random trees and gives a first glimpse at their 
fractal structure. It might also be of practical use, for instance when applied to binary 
search trees, to improve algorithmic efficiency by anticipating the structure of stored 
data. The aim of this paper is an investigation of this mass distribution, in particular 
via the derivation of scaling exponents characterizing its behavior when one side of the 
tree is much smaller than the other. Those exponents turn out to depend precisely on 
the underlying growth dynamics and also on the initial condition of the growth process. 

As long as we are interested only in global properties such as the mass distribution, 
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the growth process reduces to a model of evolution of populations which belongs to 
the class of so-called generalized Polya urn problems 0. These problems have been 
analyzed using continuous time branching processes p|-|T0[] and we believe that the mass 
distribution exponents can in principle be deduced using those methods. However, we 
find it illuminating to study the model directly using discrete time and elementary 
methods. 

We first consider growing binary trees which are characterized by a single parameter 
x given by the ratio between the probabilities that a new vertex attaches itself to an old 
vertex of valence 1 or of valence 2. For a certain value of x the model can be mapped 
onto a simple model of reinforced random walk which is exactly solvable. When starting 
from the tree consisting of a single link, one finds a uniform mass distribution, i.e. the 
probability that a fraction u of the total mass sits on the left is independent of u. More 
generally, this distribution is given by a Beta law that depends on the initial tree that 
we start the growth process from. For other, arbitrary values of the parameter x, we 
show how to calculate the mass distribution exponents exactly by scaling arguments. 
Again, these exponents depend strongly on the initial condition. 

In the second half of the paper these results are generalized to the case of growing 
trees with vertices of a uniformly bounded valence. Both for binary and the more 
general multinary trees numerical experiments have been carried out and are in excellent 
agreement with the scaling results. 

2. Growth of binary trees 

2.1. Definition of the model 

The model that we will study consists of growing binary trees from some finite initial 
binary tree T° by successive addition of links. As usual, the binary trees are rooted at 
some root edge and planar, which allows to distinguish left and right descendants at each 
vertex. One may alternatively think of the trees as drawn on top of some underlying 
infinite binary Cayley tree with a unique leaf connected to the root edge (see figure 1). 
The growing of a new link takes place at uni- or bi-valent vertices of the existing tree 
only, with two choices (left or right descendant) in the first case and a unique choice 
(the descendant not already occupied) in the second case. At each step, the choice of 
the vertex to which we attach the new link is made with a probability that depends on 
its valence i (i = 1,2) and chosen to be proportional to w± for univalent vertices and 
u>2 for bivalent ones. More precisely, if we denote by m, i = 1, 2, 3 the total number of 
z-valent vertices on the growing tree at a given step, the probability to attach the new 
link to a given i-valent vertex reads Wi/(niWi +712W2) for i = 1, 2. Note that the process 
depends only on the quantity w\/w2. Once the choice of attachment vertex is made, 
the precise choice of link to add is then unambiguous if i = 2 and uniformly distributed 
on the left and right descendants if i = 1. 

In the following, we will concentrate on global properties of the growing tree, such as 
its total numbers ni of i-valent vertices. The above growing process induces a Markovian 
evolution for ni(t) as functions of a discrete time t equal to the number of added links, 
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Fig. 1: A sample binary tree (top, thick lines) grown on top of an under- 
lying infinite binary Cayley tree (thin lines) by attaching a new link either 
at a univalent vertex (bottom left) or at a bivalent one (bottom right). 
The choice of the attachment vertex is made with probability weights u>i 
depending only on its valence i = 1 or 2. When i = 1, the new link is 
chosen among the two possible descendants with equal probabilities. 
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These rules are a particular realization of so-called generalized Polya urns with three 
types of "balls" taken out of a single urn. The problem can thus be studied by continuous 
time branching processes (see 0) but we prefer to present here some heuristic but more 
direct arguments. The rules (|2.1|) translate into a master equation for the probability 
p(ni, n 2 ,n 3 ; t) that there be rii z-valent vertices in the growing tree at time t: 



p(ni,n 2 ,n 3 ;t + 1) 



riiwx 



+ 



niwi + (n 2 - l)w 2 
(n 2 + l)w 2 



(ni - l)^i + (n 2 + l)w 2 



p(ni,n 2 - l,n 3 ;t) 



p(ni - l,n2 + 1, n 3 -l;t) . 



(2.2) 



4 



3 

This, together with the initial conditions p(rii, n 2 , 713, 0) = f] 5 n n o, determines 

p(ni,n2,ns]t) completely. 

Defining the total mass T of the evolving tree as its total number of edges (including 
the root edge), we have clearly the relation T = t + T , where T° is the mass of the 
initial tree T . Moreover, we have the standard relations for binary trees: 

T = m + n 2 + n 3 = 2ni + n 2 — 1 . (2.3) 

We may therefore write 

p(n 1} n 2 , n 3 ; t) = 5 n2jt+T o +1 _ 2ni 5 n3 , ni _i p(m; t) . (2.4) 

Hence, the problem reduces to finding a function of two variables only. 

The proportions rii(t)/T of the total mass corresponding to edges antecedent of 
z-valent vertices provide a first set of interesting quantities. As we shall see in the next 
section, these proportions tend asymptotically at large t to fixed values depending on 
W\/w2 only and not on the initial condition. Moreover, these values may be computed 
exactly via a simple mean- field argument. 

As mentioned in the introduction, another interesting quantity that characterizes 
the geometry of the growing tree is given by the repartition of the total mass between 
the left and the right of the root edge. Defining the left and right masses Tl and Tr as 
the total number of edges in the left and right descending subtrees of the root edge (with 
Tl + Tr = T — 1) , we shall be interested in the asymptotic proportions ul = Tl/(T — 1) 
and ur = Tr/(T — 1) at large t (with ur = 1 — Ul). The corresponding limiting law 
P(ul, ur) is quite subtle and cannot be obtained by mean-field arguments. As we shall 
see below, the proportions u l and ur are not peaked to fixed values but are instead 
characterized by a broad asymptotic probability distribution P(ul, ur) on [0, 1], which 
moreover depends strongly on the initial conditions. The main subject of this paper is 
an investigation of this limiting law and in particular the derivation of exponents that 
characterize it. 



2.2. Mean-field results 

The mean-field approach consists in assuming that the proportions of mass under 
study are peaked at large t around fixed values and in neglecting fluctuations around 
those. As we shall see, this gives a fully consistent result in the case of the proportions 
rii/T above, while it gives no information on the left and right proportions ul and ur. 

Multiplying both sides of eq. Q2.2Q by ni, (resp. n 2 , 77-3) and summing over all n's 
at fixed t, we get 

/ \ / \ . / n2W2 \ 
{m)t+i = (n x ) t + ) t 

niWi + n 2 w 2 
1 \ / \ . / niWl ~ n 2«^ 2 

(W2>t+1 = fa) t + t )t 2.5) 

niWt + n 2 w 2 

(n 3 )t+i = (n 3 )t + ( )t 

n\Wx + n 2 w 2 
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Here, the notation (• • -) t stands for the average over all possible growths of the initial tree 
T° by adding t links. The mean-field hypothesis allows to substitute rii — > a it within 
averages in (|2.5|) , with constants on to be determined. This substitution immediately 
yields 

a 2 w 2 otiwi - a 2 w 2 , n r , 

«i = «3 = : «2 = : (2.6) 

a\W\ + a 2 w 2 otiWi + a 2 w 2 

which fixes 

2x 2 

ai = «3 = , , ol 2 = (2.7) 

3x + a/x(8 + x) 2 + x + a/x(8 + x) 

where 

x=**. (2.8) 



a. 

i 




Fig. 2: The limiting proportions on ~ rii/t at large t as functions of the 
ratio x = 2w 2 /w\, with values (|2.7[ ) predicted by a mean field argument. 



These values are plotted in fig. 2. Note that we have the relations ol\ + a 2 + = 
2ai + a 2 = 1 in agreement with eq. (|2~3j). 

To check this mean-field result, we may use the master equation tfZ^) to estimate 
the asymptotics oip(ni,t) in (|2.4| ). Performing the substitution ( p^) into eq. (|2~^), we 
get 

p(ni; t + 1 ) = ; — : : — p(rai: t) 

FK ' ; n lWl + (t + T° + 1 - 2m)w 2 FK )! 

(t + T° + 3-2n 1 )w 2 _ 

+ (m - + (t + to + 3 - 2n t )w 2 p[Ul ' J 

At large t, we assume that the solution takes the asymptotic form 



(2.9) 



(2.10) 
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for some function f(z) and some parameter ot\ to be determined. Substituting this form 
into eq. ( |2.9|) and expanding in at large t, we get at leading order the consistency 
relation 

wi = (1 -ai)(l -2a!) i 

from which we recover the value a.\ of eq. ( |2.7|) , while the values of «2 and «3 follow 
from a 2 = 1 — 2«i and a 3 = c*i. At sub-leading order, we get a differential equation 
for f(z): 

(3 - 4ai)(/ + z/') + ai(l - ai)(l - 2c*i)/" = (2.12) 
from which we deduce that 

/(*) oc e "2-i(i-«i)(i-2«i) z2 (2.13) 



As announced, the asymptotic distribution for n±/t is peaked around the mean value a±. 
More precisely, it is Gaussian with a width of order l/y/i, hence tends to 5{n\/t — a.\) 
when t — > oo. 

It is interesting to apply the same mean-field approach to the case of the left-right 
repartition of the mass in the growing tree. The above Markov process must now be 
refined so as to keep track of the left (L) and right (R) numbers of uni- and bi-valent 
vertices n^A, i = 1,2 and A = L,R. The masses on the left and right of the root edge 
are then simply expressed as Tl = 2ni y L + ^2,l — 1 and Tr = 2n\ : R + U2,r — 1. Here 
and throughout the paper, we use the convention that n^L = 0, n 2 ,L = 1 if the left 
subtree is empty, and similarly for the right subtree. If we now assume the existence of 
limiting proportions n^A/t — > i]i,A at large t, we get the mean- field equations 

m,L = (w^L - W 2 r)2,L)/Z 

(2.14) 

m,R = {wivi,R - ^2??2,i?)/s 



where E = E^ -f- T,r and E^ = Wir]i t A + ^2??2,A- We introduce the quantities ua = 
2t\\ a + V2,A, which are nothing but the asymptotic proportions of the total mass lying 
on the left and right of the root edge. From eq. (|2.14|) , we deduce that ua = E^/E, 



henceforth ul + ur = 1 as it should. Introducing a^A = Vi,A/uA, which are the 
proportions of uni-and bi-valent vertices within the left, resp. right subtree, eqs. 02.14Q 
decouple into two sets of equations 

OCl,L = Ol 2 ,LW2/ («l,LWl + a 2y LW2) 
«2,L = (otl,LWl - a2,LW 2 )/(ai,LU)i + a 2 ,LW2) 

(2.15) 

ai,R = a2,RW 2 /(a ljR w 1 + a 2 ,RW 2 ) 
&2,r = {<*i,rWi - a 2)R w 2 )/ (a ljR W! + a 2 ,RW 2 ) 
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which are identical to the mean field equation (plf). We deduce that oq = a 



R 



cti of eq. (|2.7|) while Uj_, and ur remain undetermined. In other words, the relative 
proportions of uni- and bi-valent vertices within each subtree tend asymptotically to 
the same mean field values as in the whole tree, but the distribution P(ul,ur) is not 
peaked to particular values. We expect instead a broad distribution in the whole interval 
[0, 1]. This property will be illustrated in the next section in a particular case. 

2.3. A solvable case 

A particularly interesting case corresponds to growing the binary trees by adding 
links chosen uniformly at random among all possible available positions. This amounts 
to taking w\ = 2w2 as there are twice as many available descendants at univalent 
vertices as there are at bivalent ones. This corresponds to x = 1 in ( |2.8|) . For simplicity, 
we take it?i = 2 and u>2 = 1 so that Wim -\-ui2n2 = 2n\ +U2 directly counts the number 
of available positions for the addition of a link. This number is nothing but the mass 
T of the tree plus 1. Therefore, the evolution process no longer depends on n\ and 
individually, but instead depends on the total mass T only. 

Turning to left and right proportions, we need only keep track of the masses Tl 
and Tr of the right and left subtree (with Tl + Tr = T — 1). The induced evolution 
process consists at each step in increasing the mass by one on the left with probability 
(Tl + l)/(T + 1) and on the right with probability (Tr + 1)/(T + 1). This process is 
known as the reinforced random walk (RRW) in which a walker goes to the left (resp. 
to the right) with a probability proportional to the number of times he has already 
stepped in this direction. For a review on reinforced processes, see for instance . We 
may easily write a master equation for the probability p(TL,TR-,t) to have masses Tl 
and Tl at time t (with T = t + T°): 

p(T L , Tr; t + l) = —IL_ p (Tl - 1, T R ; t) + - + ^ - - p(T L , Tr - 1; t) (2.16) 

where T° = T° + + 1 and T® and are the initial left and right masses. Together 
with the initial condition p(Tl, Tr, 0) = S Tl T oS Tr T o , this fixes the solution to be 



p(T Ll T R ;t)= f L t+T o* t>T L +T R +i,t+T° (2.17) 



(t»)(t|) 

(t+T°\ 
\ TO ) 

In particular, when we start from an initial tree T° consisting of the root edge alone, 
we have T° = 1, T£ = = and therefore p(Tl, Tr; t) = 5T L +T R ,t/(t + 1)- Hence, all 
the values Tl = 0, 1, . . . , t are equiprobable. This is a well-known feature of the RRW. 

At large t and for arbitrary initial conditions, we may extract from eq. ( |2.17| ) the 
asymptotic law for the probability distribution of the left and right mass proportions 
ul and ur by writing P(ul, ur) = lim t 2 p(uL(t + T — 1), UrU + T — 1); t). We get 

t — >oo 

explicitly 

P(U L , Ur) = 1 yOigoT U L U R S(U L + Ur - 1) (2.18) 
L ' Ft' 
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We thus get a simple Beta law with exponents T° and the left and right masses of 
the initial tree. 

At this stage, a few remarks are in order. First, as announced, the distribution 
P(ul,ur) is supported over the whole interval [0,1]. Next, it explicitly depends on 
the initial condition. However, we note that this dependence is only through the initial 
masses and does not involve the precise shape of the initial left or right trees. In view 
of this discussion, we define the left and right exponents (3l and (3r through 

P(ul, ur) ~ u^ L when ul — > 

a (2-19) 

P(u L , u R ) ~ u R R when u R — > 

Here we have (3l = and (3r = T R , and, as the distribution is a Beta law, the left 
and right exponents characterize it completely. Note that the left exponent (3l depends 
only on the left initial condition and not on the right one and vice- versa. As we will see 
later, this property will remain true for arbitrary values of W\ and w 2 . 



3. Mass distribution exponents 

In this section, we derive the left and right exponents (3l and (3r defined in ( |2.19|) 
for the general case of arbitrary wi and 102 and for arbitrary initial conditions. 



3.1. Results 

Before we proceed to the actual derivation of the exponents, let us display and 
discuss the corresponding formula that we obtain. We have the simple expression 

/3l = — H mm {wini(T) + w 2 n 2 (T)} 

a±wi + a 2 w 2 tdT,° 

! (3-1) 

(3r = — 1 H min {wim(T) + w 2 n 2 (T)} 

aiwi + a 2 w 2 Tz>r° 



where a\ and a 2 are the limiting proportions given by eq. (|2.7|) and where the minimum 
is taken over all binary trees T containing the left (resp. right) initial subtree (resp. 
T^), with Ui(T) being the number of z-valent vertices in T. In the above formula, we 
use again the convention that whenever T is empty , we have n\(T) = and n 2 (T) = 1. 
The minimum may be explicitly evaluated and depends on the relative values of w\ and 
w 2 . Indeed, any binary tree T containing, say T®, may be obtained from by the 
successive addition of a number of links. At each step, the quantity niwi + n 2 w 2 either 
increases by w 2 > if the link is added at a univalent vertex, or is shifted by w\ —w 2 if it 
is added at a bivalent vertex. If wi > w 2 , the quantity niwi + n 2 w 2 may therefore only 
increase strictly and the minimum is attained for the initial state T = T®, with value 
w i n i l + w 2"n 2 L , where n\ L = ni(T®). When wi < w 2 , the minimum is obtained by 
saturating each bivalent vertex of into a trivalent vertex and a uni-valent one. The 
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resulting tree T has no more bivalent vertices and a number n\ L+n® L of univalent ones. 
The associated minimum now reads wi(n^ L + L ). We finally obtain the following 
explicit expression: 



0l = -1+ " V 8 ; — ' (4< L + (x + 2-\x- 2|K L ) (3.2) 

and a similar expression for the right exponent. As mentioned above, we use the con- 
vention that n\ L = and n® L = 1 if is empty. Note that the left exponent only 
depends on the left initial subtree via its numbers L of uni- and bi-valent vertices. 




X 



Fig. 3: Left mass distribution exponent (3l as a function of x = Iw^jwx 
for an initial empty left subtree, with value given by (|3]!|) with n\ L = 
and n\ L = \. At x = 1, we see the value /?£ = of section 2.3. A change 
of regime takes place at x = 2, i.e. wi = W2- 

The exponent /3c is displayed in fig. 3 for the simple case where we start from an 
empty left subtree. 

A few remarks are in order. First we note that, as a first check of our formula, 
we recover when x = 1 the result /3l — 2n\ L + n§ L - 1 = T° of section 2.3. Two 
other particular values of x may be easily checked, namely x — > and x — > oo, as the 
corresponding limiting growth processes may be easily analyzed. 

When x — > 0, we find that (3l tends to infinity unless n\ L = 0, which corresponds 
to an empty initial left subtree in which case (3l = — 1- These results may be understood 
as follows. For x — > 0, i.e. ioi >> ^2, the trees that are built are extensions of the initial 
tree by polymer-like chains (without new branching) attached to the initial leaves. If 
at least one of the left or right initial subtrees is non-empty, the number of attachment 
points for the chains is n\ L on the left and n\ R on the right, with n\ L + n\ R > 0, 
hence at large t, we have 

P{u L ,u R ) -> 5(u L q — f—Q— ) S(u R B — ~ q ) (3.3) 

n l,L + n l,R n l,L + n l,R 
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If n° L = 0, we have a delta peak at ul = 0, which is consistent with the limiting value 
(3l = —1. If L > 0, we have delta peak at a positive value of ul, consistent with a 
limiting value [3l = oo. Finally, when both initial subtrees are empty {n\ L = n® r = 0), 
the first step consist in attaching a link to the root edge with probability 1/2 on the 
left or on the right and one then grows a single chain on this side. We therefore have 
P(ul, ur) = (8(ul)8(ur — 1) + S(ul — 1)5(ur)) /2, again consistent with (3l = — 1. 




Fig. 4: A typical growth process at x = oo (from left to right). The 
process starts by first saturation all bivalent vertices of the initial tree 
T° (left) into trivalent ones, resulting in a tree S° (middle) with no more 
bivalent vertices. The growth then proceeds by a two-step elementary 
process which consist in adding the two descendant links to a leaf (circled) 
chosen uniformly among all leaves, and so on. 

Another interesting case is the limit x — > oo, i.e. w<i » w\. In this case, the 
growth process consists in first saturating all initial bivalent vertices into trivalent ones, 
resulting in a tree S° without bivalent vertex and with left and right masses S A = 
2( n i a + n 2 a) ~ 1' A = R, L (see fig. 4). Again we use the convention that n\^A = 0, 
77-2, a = 1 if the tree on side A is empty. Note that if we start from a configuration T° 
reduced to the root edge, the first two steps of the growing process consist in adding to 
the root edge the two links that descend from it. We define for consistency <S as this 
resulting tree, with Si = S R = 1. The growth then proceeds by a succession of two- 
step elementary processes that add at some existing leaf a pair of links on its left and 
right descendants. The choice of the leaf at which we attach the new pair is moreover 
uniform among all leaves. Remarkably, this two-step process with initial condition 
S° is completely equivalent to the process of section 2.3 with x = 1, where we were 
choosing single descendants uniformly. This is due to the existence of a bijection between 
"saturated" binary trees without bivalent vertices and ordinary binary trees illustrated 
in fig. 5, obtained by erasing the root edge and squeezing all pairs of descendants into 
single links. The left mass of the image of <S under the above bijection reads (#£ — 
l)/2 = n\ L + n!] ^ — 1. From the results of section 2.3, we thus get for the left exponent 
the value (3l = n\ L + n® L — 1, in agreement with ( |3.2| ) at x — > oo. 

Beyond the three solvable = 0, 1 and oo above, we can deduce a number 



of qualitative results from our general formula (|3.2|) . First, the distribution is not in 
general a Beta law. Indeed, this would require that the exponents for different initial 
conditions differ by integers, which is not the case for generic x. Second, if we start from 
the tree consisting of a single edge, the distribution can be maximal at 1^ = ^ = 1/2 
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Fig. 5: The bijection between "saturated" binary trees (top left) without 
bivalent vertices and ordinary binary trees (top right), obtained by erasing 
the root edge and squeezing all pairs of descendants into single links. The 
x = oo growth process consisting in adding both descendants of a leaf 
chosen uniformly at random (bottom left) is bijectively mapped onto the 
x = 1 growth process where single links are added uniformly at random 
(bottom right). 

for x > 1, or conversely, it can be maximal at ul = or v,r = for x < 1. Third, 
we see a change of behavior of the model at x = 2 with a change of determination for 
the exponents which suggests the existence of a competition between different effective 
growth mechanisms. Finally, the fact that the exponent 0l increases with the size of 
the initial tree is a general feature of reinforced processes. Indeed, starting from a large 
tree amounts to impose a preferential value of the mass ratio, leading to a narrower 
distribution. 

3.2. Large time behavior 

In order to estimate the, say left mass distribution exponent 0l, we are led to 
consider growth processes in which the left subtree remains finite at large time t. In 
other words, we wish to estimate the large t asymptotics of the probability p(T°, T£; t) 
to grow a tree in t steps from some initial tree T° (with left subtree T®) to a final state 
made of a fixed final left subtree T£ of finite mass T£ and an arbitrary right subtree, 
necessarily with a large mass of order t. This probability clearly decays at large t and 
we expect a power law behavior of the form 

p(T ,T£;t) ~t-T( r °' r ^) (3.4) 
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The computation of the exponent 7(T°, T£) will be a prerequisite for that of the mass 
distribution exponent (3l- 

The probability p(T°,T£;t) is the sum of the probabilities of all possible tree- 
growths of total size t leading from T° to a tree with left subtree T£. It will prove 
convenient to divide any such process into a first initial sub-process evolving up to some 
fixed time to and a later process from to to t. The value of to is kept finite but is chosen 
to be large enough so that the mean- field approximation can be applied to the (large 
enough) right subtree at each step of its evolution between to and t. We may classify 
all processes according to the intermediate configuration T° attained at time to- Note 
that T° necessarily has a left subtree intermediate between and T£ . This allows 
us to rewrite 

p{T\T£-t) = <T Q ,T Q '-M)p{T°',TL;t-t ) (3.5) 

■j-0' 

with some probabilities 7r(T°, T '; £ ) that correspond to the initial sub-process. These 
probabilities are a finite set of fixed positive numbers independent oft. 

In turn, the second part of the process may be sorted according to the induced 
left-subtree growth process. The latter corresponds to the successive addition of links 
to the left subtree until we reach T£. We will denote by Gl such a left-subtree 
growth process, characterized by the sequence of trees leading from to T£. Beside 
the precise sequence of links added, we also need to specify the time tj at which the 
j-th link is added, with j running from 1 to J = T£ — T° (where T° is the mass of 
7jP ) and with to < t± < t 2 < • ■ ■ < t j < t. Between any two consecutive such times, 
the growth process only affects the right subtree. This allows us to rewrite 

p(T°',T£;t-t )= Yl p(G L ;t 1: t 2 ,...,tj) (3.6) 

GL',tl,t2,---,tj 

where the sum extends over all left-subtree growth processes Gl with fixed initial and 
final conditions as above, supplemented by the sequence of times tj. The probability 
p{Gl'i • • • j tj) of any given such process depends only on the sequence of the numbers 
n\ = ni t L(tj) of uni- and bi-valent vertices in the left subtree at time tj (i.e., for j > 1, 
just after the addition of the j-th link). As before, we adopt the convention that n\ = 0, 
n® = 1 when is empty. The probability p(Gl; t\, . . . , tj) can be estimated as follows: 
from time tj to time tj+i — 1 (with the convention that tj+i = t), the process only takes 
place on the right, hence has a probability 

tt n liR {r)wi + n 2 ,R{r)w 2 ^ -r-r / ^ | njwi + n J 2 w 2 1 

r=4 (n{+n 1 , R (T))w 1 + (ni + n2,R(T))w 2 \ a lWl + a 2 w 2 r 

where we have used the mean-field estimate 7ij,ij(r) ~ a^r with ai as in ( |2.7|) . At time 
t = tj + i (0 < j < J — 1), the process takes places on the left, hence has a probability 
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n 



\w\ + n J 2 W2) / ((oiiWi + a-2,W2)tj+i). This results in 



/, i 1 t 2 -i _ 1 t-i _ 1 



T — t() T — t\ T — tj 

(3.8) 

where we have used the notation 

Qj = U{W1 + n * W * (3.9) 

We may now fix the left-subtree growth process Gl, i.e. keep the same sequence of link 
additions, but sum over the times tj, resulting in a probability p(Gl) such that 

p(T 0, ,T£;t-t ) = J2p( G L) (3-10) 

G L 

The corresponding nested sums over intermediate times are further approximated by 
nested integrals, with the result 

Jt <t 1 <...<tj<t V - =1 l J \ l 3 J I \ l / j=0 <ti \ l / i=0 & <?J 

(3.11) 

Again the precision of this latest approximation is governed by the choice of to, and 
can be arbitrarily accurate by picking a large enough £q. The large t behavior of p(Gl) 
therefore reads 

p{Gl) ~ where 7 = — min ^ . (3-12) 

ie{o,i,...,J} 



Introducing the notation 



for any tree T with rij(T) z-valent vertices, the minimum in ( |3.12| ) is equivalently ob- 
tained as the minimum of q(T) when T runs over all the intermediate trees of the growth 
process Gl, namely 

— min q(T) 

p(G L )~t T6G ^ . (3.14) 
Returning to p(T° , T£; t — to), we read from the summation in eq. ( |3.1U| ) that 

— min <l(T) 

p(T°',T*,t-t )~t *2' CT «2 (3.15) 

with the minimum now taken over all trees intermediate between T®' and T£. Finally, 
from eq. ( |3.5| ), we have to sum over all possible trees T ' at time to, extending the range 
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of trees in the minimum to all trees intermediate between and T£ . This results in a 
behavior p(T°,T L *;t) ~ £-T( r °. r i) with 

~/(T°,T*)= min g(T) (3.16) 
t°ctct* 

In practice, the exponent 7(T°,T^) depends only weakly on T L *. It proves useful to 
define a generic exponent 

7 (T°) = min g(T) (3.17) 
72 cr 

with no upper bound on T. In view of the discussion of section 3.1, the minimum is 
reached for some finite tree T m i n (with 7^ni n = T® if w\ > wi while, when W2 > w\, 
T m i n = S L , the tree obtained by saturating each bivalent vertex of into a trivalent 
vertex). We immediately deduce that ^(T° ,T£) takes the generic value j(T°) as soon 
as l~l contains T m [ n . Moreover, if this is not the case, we necessarily have ^(T ,^) > 
7 (T°). 

3.3. Scaling argument 

The left exponent /3l can be obtained from the exponent 7(T°) via a simple scaling 
argument. Indeed the distribution P(ul,ur) for the left/right distribution of the mass 
is obtained as 

P(u L ,u R )= lim t 2 p(u L (t + T° -l),u R (t + T° -l);t) (3.18) 

t — >oo 

where p(TL,TR)t) is the probability to have reached after t steps left and right masses 
Tl and Tr respectively. This probability implicitly depends on the initial condition T° 
(of total mass T°). We have of course Tj_, + Tr = t + T° — 1 hence we may write 

P(T L , T R ; t) = p(T L ;t) 5 Ti+THit+T o_ 1 (3.19) 

as well as 

P{u Ll u R )=P{u L ) 8{ur+ur-1) (3.20) 
where the two reduced functions are now related through 

P(u L ) = lim t p(u L (t + T°-l);t) (3.21) 

t — >oo 

The exponent (3l characterizes the behavior of P(ul) at small ul via 

P(u L ) ~ u p L L when u L -> (3.22) 



At large t and for large but finite Tl, we may according to (|3.21| ) and (|3.22|) estimate 
p{T L ;t) via 



P (T L ;t)~- t P f ~- i-f ] (3.23) 
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hence p(Ti,\t) has a power law decay t~^ L+l \ The above estimate holds in the range 
1 << Tl << t. The probability piT^t) is the sum of the probabilities p(T°,T£;t) 
over all final left subtree T£ of mass Tl- From the results of previous section, each 
p(T° ,T£;t) decays as t _ T(( r °' r i) and, as Tl is large enough, some of the T^'s will 
reach the (dominant) generic value j(T°). We deduce that 

P L + 1 = 7 (T°) (3.24) 

with 7(T°) given by ( |3.17| ). Hence, the announced result ( |3.1| ). 

3.4- Numerical checks 

To corroborate the above results, we have made a number of numerical checks. 
Those are based on random generation of large time Markov processes according to 
( j2.1|) or its left/right refinement. All simulations are made of growth processes with 
t = 4000 steps. The number iV of runs ranges from 2 x 10 5 to 10 8 according to the 
desired precision and various values of x = 2w2/w\ are explored. 
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Fig. 6: The measured distribution (see text) for n\/t (red) and n^jt (blue) 
and the corresponding mean- field predictions from ( p.7[) . 
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The first test concerns the distribution of the ratios rii/t (i = 1, 2), to be compared 
with the results of section 2.2. Figure 6 displays the measured distribution obtained 
from N = 2x 10 5 runs with an initial tree reduced to the root edge, and for x = 1/2, 1, 2 
and 4. As expected, we observe distributions peaked around precise values a.\ and a 2 
corroborating the mean field predictions (|2.7| ) . The distributions and their widths agree 
perfectly with the Gaussian shape ( |2.13| ) and its «2 counterpart. 
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Fig. 7: Measured distribution P{u£) for the initial tree T° drawn on 
the left. For x = 1, we have also indicated the exact distribution. 

Next, we have explored the left/right mass distribution for various initial conditions 
and various of x. Again the statistics is over N = 2 x 10 5 runs and we have chosen the 
values x = 1/2, 1,2 and 4. Figures 7-10 show the measured distribution P{ui) of the 
proportion ul of the total mass lying on the left side. Each point corresponds to an 
average over 10 consecutive values of the left mass, with t/10 = 400 such points. The 
initial configurations are (i) the tree reduced to the root edge (figure 7), (ii) the tree 
reduced to root edge and its left descendant (figure 8), (iii) the tree reduced to root 
edge and its two descendants (figure 9) and (iv) the tree reduced to root edge and a 
chain of length 2 on the left (figure 10). For x = 1, we have also indicated the exact 
limiting distribution as obtained by integrating ( [2.18| ) over ur. 

The very same data give access to the integrated distribution 

D{u L ) = / P{u) du (3.25) 
Jo 
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Fig. 8: Measured distribution P(u£) for the initial tree T° drawn on 
the left. 



by a simple cumulative sum. This wipes out the fluctuations of P(u), leading to 
smoother curves displayed in figure 11. In this representation, the measured x = 1 
curves are indistinguishable from their exact values deduced from ( |2.18| ). 



We now come to a check of the prediction ( |3.16| ) for the exponent 7(T°, T£) gov- 
erning the t decay of the probability p(T° ,T£;t) to go in t steps from an initial tree 
T° to a large final tree with left part T£. The predicted exponents ( |3.16| ) are plotted 
against their measured values, i.e. the limiting slopes of a log-log plot of the probability 
p(T°,T£;t) versus t. The plots 12 and 13 have been obtained from a statistics over 
N = 10 8 runs for x = 1/2 and x = 4 and starting with the initial tree T° indicated. 
The various curves correspond to various final left configurations T£ as indicated in 
medallions. For x = 1/2, all curves display the same exponent compatible with the pre- 
dicted value ( |3.17|) for the generic exponent 7(T°). Indeed, in this case, the minimum 
in ( |3.16| ) is always attained for the initial tree T° irrespectively of T£ . This is generic of 
all values x < 2. For x = 4, we see two distinct slopes according to whether T£ contains 
or not the left subtree <S° obtained by saturating the initial left subtree (see section 
3.2). If it does, we observe the generic value j{T°) ( pj.l7| ) and if not, we observe the 
larger value 7(T°,T L °) ( gJg ). 

Finally, we have made a direct measurement of the left exponent (3l from the small 
ul behavior ( |3.22| ) of P(u£). The statistics is over iV = 10 7 runs for an initial tree 
T° made of the root edge and its two descendants. The results are gathered in fig. 14 
in a log-log plot of P{u£) versus for the values x = 1/2, 1, 2,4 and 12. Each point 
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Fig. 9: Measured distribution P(ul) for the initial tree T° drawn on 
the left. 

corresponds to a left mass ranging from 1 to 100 over a total mass of 4000, hence ul 
ranges from to .025. The predicted slopes (3l are indicated by straight lines, with a 
perfect agreement in the range of large enough masses (here larger than 30) , as expected 
from the scaling argument of section 3.3. For large x, we observe a parity effect due 
to the fact that the large x preferred value n<i,L = (saturated tree) can be attained 
only for odd left masses. For x = 12, we have also represented the average between 
two consecutive left masses so as to wipe out this parity effect (empty green diamonds) . 
This seems to extend the range of validity of the scaling argument to lower it^'s. 

The very same data are used in fig. 15 to construct the log- log plot of the integrated 
distribution D{ul) ( |3.25| ) versus «l by simple cumulative sums. This wipes out the 



statistical fluctuations of the previous plot and the results corroborate unambiguously 
our predictions for (3l- 



4. Growth of multinary trees 

4-1. Definition of the model 

In this section, we generalize our results to trees with vertex valences up to k + 1 for 
some fixed integer k > 2. As before, the trees are planar and rooted and may be thought 
of as drawn on top of an underlying (fc+l)-valent Cayley tree with a unique leaf attached 
to the root edge. We start from some initial finite tree T° containing the root edge and 
add iteratively links at vertices with valence i less or equal to k, with probability weights 
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Fig. 10: Measured distribution P(ul) for the initial tree T° drawn on 
the left. 



Wi, i = 1, . . . , k. More precisely, if the tree has rij j-valent vertices (j = 1, . . . k + 1), a 

k 

given vertex of valence i is chosen with probability Wj/{ rijWj). Once the vertex is 

3=1 

chosen, the choice of link to add is uniformly distributed on all the (k + l—i) available 
descendant edges. Again, this growing process induces a Markovian evolution for the 
numbers rij(i) of z-valent vertices at time step t starting from n,(0) = the numbers 
of z-valent vertices on T . We have the evolution rules 



7li(t)tUi 



k 

I 

3=1 



m{t + l) =rn(t) 



with probability — , < ^(t + 1) = n2{t) + 1 

3=1 



m(t + l) =m(t) + l ( 41 ) 



with probability ^ J ^ + ^ = n ^ ~ 1 



„ , n i+1 (t + l) = n i+1 (t) + 1 

^nj(t)Wj { n t {t + l) =n t (t), £jtl,i,i + l 



for i = 2, . . . k. 

As before, we shall be interested in the repartition of the total mass T (number 
of edges) of the growing tree between the k descending subtrees of the root edge. If 
we denote by T m the mass of the m-th descendant subtree from the left (with T = 
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Fig. 11: Measured integrated distributions D(u) for x = 1/2 (red), x = 1 
(blue), x = 2 (cyan) and x = 4 (green) and for the four initial conditions 
indicated. 

1+X] m T m ), and w m = T m / (T— 1), we expect a broad limiting distribution P(ui, . . . , i%) 
at large T characterized by mass distribution exponents j3 m such that 



P(u!, . . . ,u k ) ~ u^" when w m -> 



(4.2) 



Mean field 

As for binary trees, the expression for the mass distribution exponents /3 m involves 
the limiting proportions oti of z-valent vertices in trees grown for a long time t. Those 
are again exactly given by a set of mean field equations 



OL\ = 1 



OilWi 



ai-iWi-i - (XiWi 
Oii = , z = 2,. 



(4.3) 



where we have defined E = Yl a i w i- These equations are easily solved into 

i=l 



Q ; 



n 



W j 



Wi L1 - E + Wj 

J= i 



(4.4) 
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3.5 4 4.5 5 5.5 
log(0 

Fig. 12: Measured exponents 7(T , T£) at a; = 1/2 from (minus) the slope 
of p(T°, T£; t) vs t in a log-log plot, for an initial tree T° as indicated and 
for various final left subtrees as shown in medallions. The black segment 
corresponds to the predicted exponent, here independent on T£. 

k 

where £ is determined by the consistency equation £ = ^ atiWi, namely 

Note that the function /(£) is strictly decreasing from /(0) = /c to /(oo) = 0, therefore 
eq. (|4.5|) has a unique real positive solution and the equation for the proportions Q4.4Q 
follow. 

4-3. A solvable case 

Here again, a particularly simple case corresponds to growing the trees by adding 
links chosen uniformly at random among all possible available positions. This amounts 
to taking Wi = (k + 1 — i) in which case ^WiUi = (k — 1)T + 1 directly counts the 
number of available positions for the addition of a link. In this case, we find E = k — 1 
and the limiting proportions oti = ( 2fc A T^^~ 1 ) / ( 2fc ^T 1 ) ■ Moreover, we can write a master 
equation for the probability p(Ti, . i) to have mass T m for the m-th descendant 
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Fig. 13: Measured exponents 7(T°,72) at x = 4 from (minus) the slope 
of p(T°, T£;t) vs t in a log-log plot, for an initial tree T° as indicated and 
for various final left subtrees as shown in medallions. The black segment 
corresponds to the two predicted exponents, depending on T£ (see text). 



subtree from the left after t steps: 



p(Tx 



,T k] t+l)= 



(k - 1)(T„ 



D + l 



^ (fc-l)(t + T°) + l 



1; 



- TO— 1 ) 



1 T 



m+1 7 



■ T fc ;t) 



(4.6) 

with initial condition p(T±, . . . , Tk] 0) = 1X71=1 ^r m ,T° • Note that this equation may 
alternatively be interpreted in terms of a generalized RRW in which a walker steps in 
one of k given directions labelled m = 1,2, . . .k with a probability proportional to a 
function f m (N m ) of the number of times N m = T m — T m he already stepped in that 
direction. Here the functions read f m (N) = N + T m + l/(k — 1). In this language, 
p(Ti, . . . , T^; t) is the probability to have stepped T m — T m times in the m-th direction 
after t steps (note that in this language T m and T m need not be integers but N m are). 
Eq.( fOD is easily solved into 



t\r(T° + T ±-) h 
p(T 1 ,...,T fc ;t) = - r 4 ^4 II 



(4.7) 
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Fig. 14: Measured exponent (3l from the slope of P(ul) vs ul in a log- log 
plot, for an initial tree T° as indicated and for various values of x. The 
straight lines correspond to the predicted exponents. 

which yields the limiting distribution P(u±, . . . , uu) = lim t k p({u m (t + T° — 1)}; t): 



P(ui, ...,u k ) 



n r (^ + ^r 



n 



t — »oo 



rjiV K-_ 

11 k - 



We can therefore read off the mass distribution exponents for this particular case: 

o _ fc-2 
Pm m fe-1 



(4.9) 



Mass distribution exponents 

In the case of arbitrary weights Wj, the mass distribution exponents may be eval- 
uated along the same lines as in Section 3, with the result 

1 k 
(3 m = -l + - min {J2 w 3 n i( T )} 



(4.10) 



3=1 



with S the solution of eq. (|4.5|) and where the minimum is taken over all trees T 
containing the m-th initial subtree 7^. Again we use the convention that whenever T 
is empty, we set ni(T) = ■ ■ ■ = nk-i(T) = and nfc(T) = 1. 
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Fig. 15: Log-log plot of the integrated distribution D(ul) vs ul from the 
same data as fig. 14. The slope are now identified with 1 +/3l, with a much 
better accuracy. 

Note that, in the solvable case of previous section, the quantity to be minimized is 
simply {k — 1)T + 1 where T is the mass of T hence the minimum is reached for 7^, 
leading to the result —1 + ((k — 1)2"^ + l)/£, in agreement with eq. ( f4.9|) as £ = k — 1. 

For illustration, let us discuss the case k = 3 in detail. It is convenient to view 
the trees T as grown out of 7^ and to follow the evolution of the quantity r(T) = 
Si=i w j n j^) that we wish to minimize. According to the evolution rules applied 
to up to three consecutive steps, r(T) has increments 

W2 when attaching one link to a univalent vertex 



(4-11) 



^3 + wi when attaching two links to a univalent vertex 

2wi when attaching three links to a univalent vertex 

W3 + Wi — W2 when attaching one link to a bivalent vertex 

2w\ — W2 when attaching two links to a bivalent vertex 

W\ — W3 when attaching one link to a trivalent vertex 

In particular, the increment is always positive whenever we attach new links to an 
existing univalent vertex. Let us introduce the ratios 



V = 



w 2 

Wi 



z = 



Wi 



(4.12) 
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Fig. 16: The four regimes (i-iv) for the determination of the mass distri- 
bution exponent (3 m in the (y — w^f-w^z = Ws/wi) plane. The exponent 
Pm is plotted here for the particular choice of 7^ indicated in the upper 
right corner. We may always write (3 m = — 1 + (ri\W\ + n 2 W2 + n^w^)/^ 
provided we choose for m the number of z-valent vertices of the appropri- 
ate tree, namely: (i) the initial tree 7^J, or (ii) a tree obtained from 7^ by 
changing its bivalent vertices into trivalent ones, or (iii) the tree obtained 
from 7^ by changing its trivalent vertices into tetravalent ones, or (iv) the 
tree obtained form 7^ by changing both its bi- and trivalent vertices into 
tetravalent ones. These trees are displayed on the right. 

For z < 1 and y < z + 1 (regime (i)), all increments above are positive. The minimum 
is therefore attained for the initial configuration T = 7^, with n^(T) = n™ the initial 
numbers of z-valent vertices in 7^. For z < 1 and y > z + 1 (regime (ii)), the quantity 
r(7~) is minimized by adding exactly one link to all bivalent vertices of 7^, thus changing 
them into trivalent ones. The resulting tree T now has ni(T) = n™ + n™, ^2(7") = 
and 713(7") = n™ + n™. For z > 1, it becomes favorable to add one link to all trivalent 
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vertices, thus changing them all into tetravalent vertices. As for bivalent vertices, either 
y < 2 (and in particular y < z + 1) and we leave these vertices unchanged (regime (iii)), 
resulting in a tree with ni(T) = n™ + n™, ^(T) = n™ and n^{T) = 0, or y > 2 and 
we must change these bivalent vertices into tetravalent ones (regime (iv)), resulting in 
a tree with ni(T) = n™ + 2n™ + n™, n2(T) = n^{T) = 0. This allows us to write 

(i ) /jm = "r-i+"^+^3 _ 1 for z<hy<z + 1 

(ii) ft. = jPl±^l±M±^l _ ! for z<1 , y>z+1 

E (4 13) 

(iii) Pm = — ^ — - 1 for z > 1, y < 2 

(iv) (3 m = ^ 1 - 1 for z > 1, y > 2 

These formulae also hold whenever 7^ is empty by setting n™ = n™ = and n™ = 1. 
These results are summarized in figure 16 where the different regimes are displayed. 
We have also performed a number of numerical checks identical to those presented in 
section 3 and which fully corroborate the above results. 



5. Conclusion 

In this paper, we have derived the exact expression of exponents characterizing the 
mass distribution of trees growing locally by addition of links. As opposed to generic 
random trees whose statistics is governed by local Boltzmann weights, the mass of large 
growing trees is not in general concentrated in a single side of the root but can be 
distributed in all subtrees. The mass can be preferentially concentrated in one subtree 
(negative mass distribution exponent), it can be preferentially equally distributed over 
all subtrees (positive exponent) or it can be characterized by a uniform distribution 
(vanishing exponent). In practice, the precise value of the mass distribution exponent 
depends both on the growth process parameters Wi and on the initial condition T . In 
particular, in the case of binary trees grown from a single edge, the exponent Pl increases 
with x up to the transition point x = 2 and then decreases down to at x —>■ oo. This 
shows that a monotonous variation of the relative strength that we attach to the two 
types of vertices (univalent or bivalent) does not induce a monotonous variation of the 
characteristics of the tree. 

We have observed that, except for special values of x where the model is exactly 
solvable, the mass distribution is not in general a Beta law and we may wonder whether 
an explicit analytic expression could be found. In this respect, a promising direction of 
investigation is provided by the approach of the problem via continuous time branching 
processes flSHlQH- Indeed, within this framework, it is possible to relate the mass distri- 
bution to other asymptotic large time distributions which are themeselves determined 
by coupled integral equations. Solving these equations analytically would be a major 
advance in this field. 
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For more general tree valences, an interesting outcome of our study is the exis- 
tence of different regimes for the value of the mass distribution exponents. This seems 
to indicate the existence of a competition between various coexisting effective growth 
mechanisms. Whether the transitions between these regimes are characterized or not 
by the emergence of some order parameter remains to be understood. 
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